function cfdInitialReadSteadyRestartFile
%--------------------------------------------------------------------------
%  Written by LiuHaHa @ NWPU, 2022 01 
%  Contact me at: liuzhikan@mail.nwpu.edu.cn
%==========================================================================
% Routine Description:
%   This function initialize the flow field using steady restart data file 
%--------------------------------------------------------------------------
global cfdFlow slvpara;
phyTemp = cfdFlow.Phy;

filename=[slvpara.meshFile,'/Restart/','steadyRestart.dat'];
fid=fopen(filename,'r');

tempA = mod(size(phyTemp,1),3);
tempB = round((size(phyTemp,1) - tempA)/3);
for numVar=1:4
    for item=1:3:tempB*3    
        [phyTemp(item:item+2,numVar)] = fscanf(fid,'%f %f %f ', [3 1]);   
    end
    tempC = size(phyTemp,1) - tempB*3;
    if (tempC == 1)
        [phyTemp(tempB*3+1:tempB*3+tempC,numVar)] = fscanf(fid,'%f ', [tempC 1]); 
    elseif (tempC == 2)
        [phyTemp(tempB*3+1:tempB*3+tempC,numVar)] = fscanf(fid,'%f %f ', [tempC 1]); 
    end
end

[cfdFlow.vel(:)] = fscanf(fid,'%f %f', [2 1]);

fclose(fid);    

cfdFlow.Phy = phyTemp;
















